Solitons in a medium with linear dissipation and localized gain 
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We present a variety of dissipative solitons and breathing modes in a medium with localized 
gain and homogeneous linear dissipation. The system posses a number of unusual properties, like 
exponentially localized modes in both focusing and defocusing media, existence of modes in focusing 
media at negative propagation constant values, simultaneous existence of stable symmetric and 
anti-symmetric localized modes when the gain landscape possesses two local maxima, as well as the 
existence of stable breathing solutions. 

PACS numbers: 42.65.Tg, 42.65.Sf 



Dissipative solitons appear in a variety of forms and 
settings. They form in laser systems with saturable 
gain and absorption [l[ , in homogeneous materials where 
light evolution can be described by the cubic-quintic 
Ginzburg-Landau (GL) equation 0, and in systems with 
spatially localized gain and uniform nonlinear losses 0- 
^ . A characteristic feature of all these systems is that the 
existence of dissipative solitons there is mediated by the 
competition of gain and losses having different physical 
origin. Thus, in the cubic-quintic GL equation solitons 
exist because of exact balance of the linear losses, cubic 
gain, and quintic losses. Likewise, solitons can also form 
because cubic losses compensate the spatially localized 
gain. This has been established in one-dimensional set- 
tings in the presence of a linear localized potential Q or 
lattice , in two-dimensional problems in planar waveg- 
uide arrays Q , or lattices @ , or even in two-dimensional 
systems without additional refractive index modulations 



It is therefore important to elucidate simpler physical 
settings, where only the linear gain and losses result in 
formation of dissipative solitons. This is the goal of this 
Letter. More specifically, we report on the existence of 
stable dissipative solitons supported by a localized gain in 
a homogeneous Kerr medium where only the linear losses 
are present far from the center of gain landscape. This 
simple system possesses a number of interesting prop- 
erties among which is the existence of localized modes 
even for negative propagation constant values, simulta- 
neous stability of symmetric and antisymmetric modes 
in two- humped gain landscapes, as well as the existence 
of breathing modes. 

We consider the nonlinear Schrodingcr (NLS) equation 
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(1) 



where a = 1 (cr = — 1) corresponds to defocusing (focus- 
ing) nonlinearity, j{x) describes a localized gain concen- 
trated in the domain (or domains) having a characteristic 
width d. Without loss of generality we rescale the back- 



ground dissipation to the value: 70 = 1, and consider the 
situation where max7(x) > 1. 

We look for spatially localized solutions obeying 
lim|3.|^oo |(7(z, x)| — and start with the stationary 
modes q{z,x) = e^'"^w{x) where w{x) = e^^^'^'^u[x) with 
u{x) > and 9(x) being real functions. Then, introduc- 



ing the current density i{x) 



one obtains from ^ 



the system describing stationary modes 



bu 



0, j^ + [l-^{x)]u^ (2) 



In order to establish the domain of possible variation of 
the propagation constant b, we now address the limit 
X ^ 00 focusing on the exponentially decaying solutions, 
i.e. assuming u ^ Ce^^^, where C and ^ are positive 
constants. From the second of Eqs. ^ it readily follows 



that in this limit j ^ C 



/2^. This suggests that in 



the first of equations ^ the term au^ can be neglected 
in the limit x — > cx), while substituting the asymptotic 
values of u and j in the terms which are left we obtain 
b = fi^ — l/(4/i^). It follows from this formula that the 
propagation constant for localized solutions can be either 
positive or negative even if the nonlinearity is focusing 
(contrary to what happens in uniform conservative sys- 
tems or in dissipative systems 0, 0, S, 01)- If ^ ^ +00 
then /i — > -l-oo while if 6 —00 then ^ +0. Respec- 
tively, there must exist localized modes with /i > even 
at 6 = while the spatial delocalization occurs in the 
limit b — !■ —00. 

The condition of balance between the dissipation and 
gain (for 70 = 1) has the form U ~ J j{x)\u\'^dx where 
U — J |up is the energy flow. From the balance be- 
tween the nonlinearity and diffraction we obtain 



bU=- 



'dx 



'\^dx 



(3) 



From this formula one readily concludes that for zero 
propagation constant, 5 = 0, the mode should decay with 



2-1/2 



I.e. 



^ 0, and hence ^ can be satisfied 



2 



only in the focusing medium (i.e. at cr = —1). Thus 
the linear limit of the problem (i.e. U ^ 0) corresponds 
to a localized mode with b = 6*^°) < (indeed, in this 
case J \w\'^dx/U < maxjwl^ — >• and can be neglected 
in ©). In the defocusing medium exponentially local- 
ized solutions exist only if & < 6^"^ while in the focusing 
medium dissipative solitons can exist with both positive 
and negative propagation constants, i.e. at 6 > 

We start with the case of the gain localized in one 
spatial domain, setting 7(0;) = ae~^ , (a > 1 char- 
acterizes the gain strength) . By fixing the width of the 
gain landscape d and changing a we obtain the depen- 
dencies b{a) and U{a), which are depicted in Fig. [1] for 
d = 1. These dependencies bifurcate from the linear limit 
corresponding to U = and to b'"^ w —0.61. The very 
fact of existence of localized linear modes with b — 6*^°' 
in inhomogeneous gain landscapes is known as the gain- 
guiding effect (see e.g. 0). In the linear case the guiding 



occurs for a particular value of a 



7.(0) 



1.98. Focusing 



nonlinearity diminishes gain at which localization occurs, 
while in a defocusing medium a higher gain is required 
for localization (Fig. [IJ. In Fig. [1] we also show an ex- 
ample of a nonlinear mode. Increase of the gain results 
in decrease of U in the focusing medium, that is also ac- 
companied by expansion of the soliton. In the defocusing 
medium the growth of a results in increase of the energy 
flow and progressive expansion of the soliton outwards 
the gain domain. Solitons are exponentially localized in 
the both focusing and defocusing media, as predicted by 
the analysis of ([2]) which showed that the current results 
in effective renormalization of the propagation constant. 
Physical understanding of the exponential localization in 
the defocusing medium, comes also form the observation 
that the current results in the effective self-focusing of 
the beam, what stems form the negative sign in front of 

Linear stability analysis (confirmed by the direct prop- 
agation) indicates that the solitons in the focusing 
medium are unstable [see Fig. Elja)], while in defocus- 
ing medium they are attractors and can be excited with 
a variety of inputs ranging from noisy to localized pat- 
terns. This difference in the stability can be understood 
from Eqs. The solution represents a fiow outwards 
the "source" (i.e. the gain domain). The defocusing 
nonlinearity enhances the outflow from the high inten- 
sity region thus contributing to the stability, while the 
focusing medium enhances the field concentration in the 
gain domain, what stimulates further growth of the peak 
amplitude. 

Next we study the case of the gain landscape with two 
maxima at a; = ±xo considering "f{x) = ae^(^~^°'' Z'' + 
fj_^-{x+xo) Id ^ .^^rijgj-g a > 1. Here one also can obtain 
modes both for focusing and defocusing media. In Fig. [3] 
we present symmetric (even) and antisymmetric (odd) 
modes. Both types of solutions bifurcate from the linear 
limit. The values of parameters a and b at which bif urea- 
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FIG. 1: Propagation constant (a) and energy flow (b) vs a 
for solitons supported by single gain channel. Red and black 
curves correspond to focusing and defocusing media. Dashed 
line indicates a value corresponding to the linear limit. Profile 
of soliton (c) and current density (d) in defocusing medium 
at a = 2.5 corresponding to the circles in (a) and (b). 
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FIG. 2: (a) Instability increment 5r vs a for solitons in a sin- 
gle gain channel in focusing medium. Circles indicate points 
beyond which solitons exist only in defocusing medium where 
they are stable, (b) Sr{a) for symmetric "s" and antisymmet- 
ric "a" solitons in two gain channels with d = 1 in defocusing 
medium. 
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tions take place are very close (6 
for symmetric, and fei^' w —0.60, ai*^' w 1.98 for anti- 
symmetric solitons). Closely to the linear limit and in 
the focusing medium the U{a) curves for the symmetric 
and antisymmetric solitons are almost indistinguishable, 
while in the defocusing medium their energy fiows may 
differ considerably [Figs. [3] (a) and (b)]. The examples of 
modes are shown in Figs. [3] (c) and (d). 

Linear stability analysis reveals multiple alternating 
stability and instability domains (in a) for both symmet- 
ric and antisymmetric modes [Fig. [2] (b)]. A remarkable 
fact is that stable symmetric and antisymmetric modes 
co-exist at the same parameters of the system. Moreover, 
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FIG. 3: Propagation constant (a) and energy flow (b) vs a 
for solitons supported by the two gain channels at so = 6. 
Solid and dashed lines correspond to the symmetric and an- 
tisymmetric modes. Red and black lines correspond to the 
focusing and defocusing media. The circles correspond to the 
solitons shown in (c),(d). (c) Unstable antisymmetric soliton 
in the focusing medium at a = 1.8. (d) The stable symmetric 
soliton in defocusing medium at a = 2.6. 

for the same value of a both symmetric and antisymmet- 
ric modes in the defocusing medium can be attractors. 

We found that for certain values of a even input beams 
evolve into the symmetric stationary mode while odd in- 
put beams evolve into the antisymmetric mode, see pan- 
els (a) and (b) of Fig. |4l When one of the modes is 
unstable it usually evolves into a stable mode or into pul- 
sating mode. An example of such behavior is shown in 
Fig. 131(c) where perturbed unstable antisymmetric mode 
transforms into the breather. (We show only the initial 
stage of the evolution; simulations were performed up to 
z = 2 • 10^). Dynamics shown in Fig. |4l^c) is also illus- 
trated in Fig. [5] where a corresponding trajectory in the 
phase space and dependence U (z) are plotted. 

Finally we explore the possibility of switching between 
stable solutions by changing the control parameter a. An 
example is shown in Fig. 2] (d) where the input beam 
corresponds to the stable symmetric mode at a = 2.186. 
At z — 100 the strength of the gain a is switched from 
2.186 to 2.5. Since the symmetric mode at a = 2.5 is 
unstable, the switching results in the beam evolution into 
the stable antisymmetric mode corresponding to a = 2.5. 

To conclude, we found a variety of unusual properties 
of the solitons described by model ([T]) , such as existence 
of exponentially localized modes with negative propa- 
gation constant values in both focusing and defocusing 
media, simultaneous stability of symmetric and antisym- 




FIG. 4: Excitation of stable symmetric (a) and antisymmet- 
ric (b) solitons by symmetric and antisymmetric input beams 
in a system with two gain channels for a = 2.391. (c) Ex- 
citation of a breather at a = 2.61 starting with the unstable 
antisymmetric soliton. (d) Switching between stable symmet- 
ric and antisymmetric solitons stimulated by abrupt change 
of a at 2: = 100. 




FIG. 5: The trajectory (panel a) and dependence U{z) (panel 
b) corresponding to the dynamics shown in Fig. [IJc). The 
bold line shows the limit cycle. 



metric solitons, and the existence of persistent breathers. 
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